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Abstract 

The quark-meson coupling (QMC) model, which has been successfully used to describe 
the properties of both infinite nuclear matter and finite nuclei, is applied to a systematic 
study of A, S and H hypernuclei. Assumptions made in the present study are, (i) the 
(self-consistent) exchanged scalar, and vector, mesons couple only to the u and d quarks, 
and (ii) an SU(6) valence quark model for the bound nucleons and hyperon. The model 
automatically leads to a very weak spin-orbit interaction for the A in a hypernucleus. Ef- 
fects of the Pauli blocking at the quark level, particularly in the open, coupled, S A — AA 
channel (strong conversion), is also taken into account in a phenomenological way. 
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1 Introduction 



In earlier work [|I], 0, |^ we addressed the question of whether quarks play an important role in 
finite nuclei. This involved quantitative investigations of the properties of closed-shell nuclei 
from ^^O to ^ospb |l|, as well as the effective mass of the p meson formed in light nuclei 
These calculations were performed within the quark-meson coupling (QMC) model, originally 
suggested by Guichon [|], and with an extended version Q], which treats also the u and p 
meson mass variations self-consistently in the medium. 

In the QMC model, the interactions between hadrons are mediated by the exchange of 
scalar (cr) and vector {u and p) mesons self-consistently coupled to the quarks within those 
hadrons. Within the model it has proven possible to successfully describe the properties of 
infinite nuclear matter ^ and also finite nuclei ^ ^. Blunden and Miller and Jin 
and Jennings have made similar studies based on the QMC model, and some phenomeno- 
logical extensions. One of the most attractive features of the QMC model is that it does not 
involve much in the way of additional complications to Quantum Hadrodynamics (QHD) ||^. 
Furthermore, it produces a reasonable value for the nuclear incompressibility IQ-IQI- 

Here we apply the QMC model [|I], ^ to a systematic study of the properties of A, S and S 
hypernuclei. Some of the initial results for A hypernuclei were reported already 0. One of the 
purposes of the present article is to present the spin-orbit potentials for the hyperon calculated 
self-consistently with an explicit quark structure for the bound hyperon in the QMC model. 
These spin-orbit potentials generally contain anomalous contributions from the quarks due to 
the finite size of the hyperon. 

Within the Born-Oppenheimer approximation, one can derive equations of motion for a 
hypernucleus in the QMC model in the same way as has been done for normal nuclei 
H. Such an approach may provide us with important information about the hyperon-nucleon 
interaction, the deep nuclear interior and a possible manifestation of the quark degrees of 
freedom via the Pauli principle at the quark level |]10|-|^- Below, we will briefiy review the 
present situation regarding A hypernuclei first, because they have been well studied, both 
experimentally and theoretically. 

As an example, the very weak spin-orbit interaction for A hypernuclei, which had been 
phenomenologically suggested by Bouyssy and Hiifner |jlO[, was first explained by Brockman 
and Weise in a relativistic Hartree model, and directly confirmed later by experiment [|I^ . 
However, a very strong SU(3) breaking effect was required to achieve this small spin-orbit force. 
An explanation in terms of quark and gluon dynamics was made by Pirner |13||. Alternatively, 



Noble showed that this small spin-orbit force could be also realized without any large 
breaking of SU(3) symmetry, if an tuAA tensor coupling was introduced, analogous to the 
anomalous magnetic moment of the A. However, Dover and Gal [jl5| questioned whether the 
tensor coupling of the uj meson to the A could be related to the anomalous magnetic moment, 
because the spin of A is entirely carried by the s quark in a naive SU(6) valence quark model, and 
this s quark couples exclusively to the meson according to the OZI rule. Later, Jennings 
pointed out that within Dirac phenomenology the tensor coupling of the uj meson to the A 
could be introduced in such a way as to guarantee that the direct uj coupling to the spin of the 
A was zero - as one would expect from a simple quark model. The resulting spin-orbit force 
agrees with the result obtained in the present model |p, on the basis of an explicit treatment 
of the quark structure of the A moving in vector and scalar fields that vary in space. 

As a second example, it has also been discovered that there is an overbinding problem in 
the light A hypernuclei, and the existence of a repulsive core or the necessity of a repulsive three- 
body force have been suggested to overcome the problem ||T^ . The origin of this overbinding. 
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which could not be explained easily in terms of traditional nuclear physics, was ascribed to 
the Pauli principle at the quark level by Hungerford and Biedenharn [0. Investigations of 
this repulsive core in the A-nucleus system have been made by Takeuchi and Shimizu, and 
others |T9|, ^ , based on a nonrelativistic quark model. 

In addition to the investigations of A hypernuclei, many studies of S and S hypernuclei 
have been also made 111,11-11,10 



|. In particular, experiments on E hypernuclei per- 
formed at CERN , Brookhaven and KEK , confirmed the existence of S hypernuclei. 
The former two experiments were done using K~ beams in so called "in-flight kinematics", 
while the latter was done using stopped kaons. These experiments revealed an interesting fea- 
ture of S hypernuclei, namely, the widths of the S in a hypernucleus seems to be narrower 
than the naive expectations due to the strong EiV — AA^ conversion (channel coupling). As for 
the S hypernuclei, although not so many experiments on the double strangeness hypernuclei 
(S = —2) have been performed so far |BD]-|BB], there is a renewed interest in connection with 
the H dibaryon suggested by Jaffe ||6^, based on quark degrees of freedom - using the MIT 
bag model. Thus, we expect that more experiments on the S and S hypernuclei, with higher 
statistics and precision, will be performed in the near future. 

In the light of these investigations, it now seems appropriate to investigate the (heavier) 
hypernuclear systems. A, S and H hypernuclei quantitatively, using a microscopic model based 
on quark degrees of freedom. For this purpose, the QMC model (which is built explicitly on 
quark degrees of freedom) seems ideally suited, because it has already been shown to describe 
the properties of finite nuclei quantitatively. The present investigation is one of the natural 
extensions of the study of quark degrees of freedom in finite nuclei. 

The organization of the paper is as follows. In Section 2, the relativistic formulation of the 
hypernuclear system in QMC will be explained. A mean-field Lagrangian density and equations 
of motion will be derived in Section 2.1. Hyperons in infinite nuclear matter will be discussed in 
Section 2.2, while (finite) hypernuclei will be treated in Section 2.3. In Section 3, the spin-orbit 
potential for the hyperon in the QMC model will be discussed. Some results for hypernuclei 
in the QMC model alone will be given in Section 4. In Section 5, modifications necessary 
for a more realistic calculation in addition to the present version of the QMC model will be 
discussed. Effects of the Pauli blocking at the quark level in the AA^ — UN (and "EN — AA) 
coupled channel, will be considered in specific ways at the hadronic level in Sections 5.1, and 
5.2, respectively. The results with including these effects discussed in Sections 5.1 and 5.2 will 
be given in Section 5.3. Finally, Section 6 will be devoted for summary and discussion. 



2 Hypernuclei in the QMC model 

In this Section, we will derive mean-field equations of motion for a hypernucleus, as well as 
consider a hyperon in nuclear matter. 

2.1 Mean-field equations of motion 

Using the Born-Oppenheimer approximation one can derive mean- field equations of motion 
for a hypernucleus in which the quasi-particles moving in single-particle orbits are three-quark 
clusters with the quantum numbers of a hyperon or a nucleon. One can then construct a 
relativistic Lagrangian density at the hadronic level 0, |^, ^ , similar to that obtained in QHD |^ , 
which produces the same equations of motion when expanded to the same order in velocity, v: 

rHY _ r , rY 

'-QMC — '-QMC -r f-QMCi 
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C-QMC = i^N^r) 



3 I /I 1 



+ limf))' + mib{fr] + ^(vA(f))^ 



y=A,s,H 

where ipNi'f^) ('?/'y(r*)) and 6(r*) are respectively the nucleon (hyperon) and the p meson (the 
time component in the third direction of isospin) fields, while rrifj, and rrip are the masses 
of the a, uj and p mesons. (7^ and gp are the w-N and p-N coupling constants which are related 
to the corresponding (u,d)-quark-ci;, g^, and (u,d)-quark-p, (7^, coupling constants as g^ = 3g^ 
and gp = g1, §,§. 

In an approximation where the a, uj and p mesons couple only to the u and d quarks 
(ideal mixing of the u and (p mesons, and the OZI rule are assumed), the coupling constants in 
the hyperon sector are obtained as g^ = {no/3)guj, and gj = gp = g"^, with uq being the total 
number of valence u and d quarks in the hyperon Y. and Qy are the third component of the 
hyperon isospin operator and its electric charge in units of the proton charge, e, respectively. 
The field dependent cx-N and cx-Y coupling strengths predicted by the QMC model, g^icr) and 
g^{cr), related to the Lagrangian density Eq. (0) at the hadronic level, are defined by: 

M*(a) = M^-g,{a)a{r), (2) 
M^(a) = My-g^ia)air), (3) 

where Mn {My ) is the free nucleon (hyperon) mass. Note that the dependence of these coupling 
strengths on the applied scalar field must be calculated self-consistently within the quark model. 
Hence, unlike QHD, even though g^{<y)/ga{o') may be 2/3 in free space (cr = 0)[] , this will not 
necessarily be the case in nuclear matter. More explicit expressions for g^{(y) and 5'o-(cr) will 
be given later. From the Lagrangian density Eq. (|T|), one gets a set of equations of motion for 
the hypernuclear system: 

[i^-d- Ml,{a) - ( gM^ + gp^Kv) + -(1 + T^)A{f) )7o]^^(r) = 0, (4) 
[i^-d- M^{a) - ( glcuir) + gpljb{f) + eQyA{r) )7o]^y (r) = 0, (5) 
{-Vl + ml)a{r) = -[ > ,(r-) - [ > r(r), 

= g.CN{<y)p,{T) + glCy{a)p^,{r), (6) 

{-Vl + ml)uj{f) = g^pBif) + g^pli^, (7) 

i-Vl + ml)b{r) = ^p,{r) + g^^ pl{f^, (8) 

{-yl)A{T) = ep,{r) + eQyplir), (9) 

where, Ps(r) (pf (^)), Pb(^) {Pb{^): P3(^) and Pp(r) are the scalar, baryon, third component of 
isovector, and proton densities at the position r in the hypernucleus p|, H]. On the right hand 



side of Eq. (g), a new, and characteristic feature of QMC beyond QHD g, |3^, ^ appears 



"'^ Strictly, this is true only when the bag radii of nucleon and hyperon are exactly the same in the present 
model. See Eq. ([l^), below. 
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namely, 
9l{^ = 



gaCN^cr) and 



(9Af*((T) 



g^Cvifr), where = g^{a = 0) and g] 



0). The effective mass for the hyperon Y is defined by 

My(cr) = — h ■^TT{nY) B, 



with the MIT bag model quantities |1|, |^ 



Cy(a) 
9^ 



Syia)/SYiO), 

noglSYiO) = ^fi'a5'y(0)/S'iv(0) = ^g^TY/N 



[a 



^x^ + {RYml,{a)y 
mu4- gl(y{r), 



x^ + {RYrris 



such Ry that satisfies 



dRy 



0. 



The density dependent a-Y coupling constant is defined by 



dM*{a) 
da 



bag 



- d 

i^uAy)^uAy) = -n^giSYio) = -— 



(10) 



(11) 

(12) 
(13) 

(14) 
(15) 

(16) 



(17) 



For the nucleon, replace the suffix F by in Eqs. (|T0|) - ( [T7| ) and set no = 3. Here, zy, B, x, Xg 
and mu^d,s are the parameters for the sum of the cm. and gluon fluctuation effects, bag pressure, 
lowest eigenvalues for the (u, d) and s quarks, respectively, and the corresponding current quark 
masses with m„ = rrid. zn and B (zy) are fixed by fitting the nucleon (hyperon) mass in free 
space. The bag parameters calculated and chosen for the present study are, Rn = 0.8 fm (in 
free space), B = (170 MeV)^ m„ = = 5 MeV and m, = 250 MeV. Note that due to the 
equilibrium condition Eq. (^), which is applied also to the free nucleonic vacuum, the bag 
radius for each hyperon in free space is slightly different. The parameters zj^j^y and the bag 
radii, Rn,Yi obtained by fitting to their physical masses in free space are listed in Table j^. 



Table 1: The bag parameters, -Zat^y, the bag radii in free space, Rn.y and the physical baryon 
masses fitted in free space. They are obtained with the bag pressure, B = (170 MeV)^, current 
quark masses, rriu = rrid = 5 MeV and = 250 MeV. 





mass (MeV) 


zn,y 


Rn,y 


(fm) 


N 


939.0 


3.295 


O.J 


^00 


A 


1115.7 


3.131 


O.J 


^06 


S 


1193.1 


2.810 


O.J 


^27 




1318.1 


2.860 


0.5 


^20 



The parameters associated with the u and d quarks are those found in our previous 
investigations [^. The value for the mass of the s quark was chosen to be 250 MeV. We note 
that the final results are insensitive to this parameter. The value for Ty/n in Eq- (0) turned 
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out to be almost unity for all hyperons even though Rn 7^ Ry (see also Eqs. (|TT])), so we can 
use Ty/N = 1 in practice ||^. Although the bag radii in free space for the nucleon and hyperons 
are slightly different {R^ 7^ Ry), they will be generically denoted by Rb = 0.8 fm in Figs. |I] 
and 

At the hadronic level, the entire information on the quark dynamics is condensed into 
the effective coupling CN,yicr) of Eq. (H), or Eq. ([T2|) . Furthermore, when CN,yio') = 1, which 



corresponds to a structureless nucleon or hyperon, the equations of motion given by Eqs. 
can be identified with those derived from QHD [^, 



42| , except for the terms arising from 



the tensor coupling and the non-linear scalar field interaction introduced beyond naive QHD. 



2.2 Nuclear matter limit 

Here, we consider a hyperon in nuclear matter. In this limit the meson fields become constants, 
and we denote the mean- value of the a field as W. The self-consistency condition for the a field, 
a, and effective nucleon mass, M^, is given by ||, || 



a = ^C'iv(o-) 777-7:7 / dk9{kF - k) 
mi ivr p J 



(2^)' J \lMf{a) + P 



'TV 

where = {3g^SN{0)), kp is the Fermi momentum, and CAr(a) is now the constant value 
of Cat in the scalar field. Note that Mj^(a), in Eq. (|TH|), must be calculated self-consistently 
by the MIT bag model, through Eqs. ([10|) - ([T7|) . This self-consistency equation for W is the 
same as that in QHD, except that in the latter model one has Cat (a) = 1 |^. By using the 
obtained mean field value, a, the corresponding quantity for the hyperon Y, Cy(a), can be 
also calculated using Eqs. ( |lTl) and (|T2l), where the effect of a single hyperon on the mean field 
value, a, in infinite nuclear matter can be neglected. The calculated values of Cj^y{a) for the 
nucleon and hyperons are displayed in Fig. |^. It has been found that the function Cjqy{a) can 
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Figure 1: C{a) for nucleon and hyperons in nuclear matter. 

be parametrized as a linear form in the a field, gad", for practical calculations [|l], H, H by 

CN,y{^) = 1 - a-NX X {gao), (19) 
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which can be also seen from Fig. |T], where g^ja = (3(y'^S'Ar(0)) a in MeV. The values of oat y for 
nucleon and hyperons are given in Table |^. This parametrization works very well up to about 



Table 2: The slope parameters, and ay. 





xlO-^ MeV-i 




8.8 


OA 


9.3 


as 


9.5 


Oh 


9.4 



three times normal nuclear density, pb — Bpo, with po — 0.15 fm"'^. For the field strength, gc^a, 
versus baryon density, see Ref. [^. 

Next, we define the variations of effective masses for the baryon j {j = N, A, S, S): 

6M* = Mj - M;, (20) 

which should be calculated self-consistently through the corresponding Eqs. ([T0| ) -([18| ) . Calcu- 
lated values for the variation of effective masses against baryon density are shown in Fig. ^ 
For the effective mass, M*, one can write down explicit expressions using the parametrization 



350 




0.0 0.5 1.0 1.5 2.0 2.5 



Figure 2: Variation of the effective masses for nucleon and hyperons in nuclear matter, where 
normal nuclear matter density, po, is 0.15 fm~^. 

for C^,y(a) of Eq. (|19|): 



m: 



N 



Mn - ga{(y)(y ^Mn - g^ 



MY-gl{a)a^MY-'jgo 



1 - '^i.9<rcr] 



cr. 



(21) 
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Because the parameters aj^y are almost the same in magnitudes (See Table we find a simple 
scaling relation for the variation of effective masses [H: 



6MI ^ dM* ^ 2 5M; _ 1 
6M% ~ SM^ ~ 3' 6M^ ~ 3' 



(22) 



This scaling relation is well demonstrated in Fig. 0. It is worth noting that the ratio of the 
effective mass to the bare mass, M^/Mj, does not decrease linearly as the baryon density 
increases in the QMC model, which is different from the relation adopted by Kuwabara and 
Hatsuda W^, based on a model of the QHD type. 



2.3 Hypernuclei 

The variation of the Lagrangian density Eq. results in the following equations for static, 
spherically symmetric nuclei plus one hyperon configuration (Hereafter, we will call this con- 
figuration simply a hypernucleus) : 

^2^ir) + ^^air)-mlair) = -g.C^{a{r))ps{r) ~ glCy{a{r))p': {r), (23) 
d? 2d 

^^l'^) + - ^l^i^) = -9cvpBir) - g^plir), (24) 

;b{r) + --^b{r) - mjb{r) = -^psir) - glt^plir), (25) 



where 



^ 1 / \ ^ ^ T / \ 2t/\ Q P / \ Y I Y / 



^^2^W + -^^W = -ePp{r)-eQypUr), (26) 

(27) 



(28) 
(29) 
(30) 
(31) 
(32) 

(33) 



Psir) 


= E^«' 

a 


:r)(|G.(r)p-|F.(r)r), 


Pjir) 


= dp{r){ 


\G,{r)\'~\F,{r)\'), 




occ 




pB{r) 


= E^°' 

a 


[r){\G^{r)\' + \F^{r)\% 


Plir) 


= dp{r){ 


\G,{r)\' + \F,{r)n 




occ 




Psir) 


= J2da^ 

a 




Ppij) 


occ 

= E^"' 
a 


[r){ta + l){\G^{r)\'+\F^{r)\') 



with da{r) = {2ja + 1)/A7rr^, dp{r) = l/Airr^ and 

d K, 1 

:rGa{r) + -Ga{r) - [e^ - g^r) - tagpb{r) - (t„ + -)eA{r) + Mn 
dr r 2 



d K 1 

-^F^{r)--F^{r) + [e^-g^uj{r)-t^gpb{r) - (t„ + -)eA(r) - Mjy 



(7,(a(r))a(r)]F,(r) = 0, (34) 
(t„ + ^)eA(r) - Mjv 
+ g,{a{r))a{r)]G^{r) = Q, (35) 
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eQyv4(r) + My 

gl{a{r))a{r)\ Fp{r) = 0, (36) 
eQyA(r) - My 



0. 



(37) 



Here Ga,i3{r)/r and Fa^i3{r)/r are respectively the radial part of the upper and lower components 
of the solution to the Dirac equation for the nucleon (hyperon) 



i^N,Y{r) 



2[G,,^(r)/r]<l> 
-[F,,^(r)/r]$ 



Km 
Km I 



(38) 



where ^ is a two-component isospinor and is a spin spherical harmonic for the nucleon 
(hyperon) (a (/?) labelling the quantum numbers of nucleon (hyperon) and e^^p being the 
energies). Then, the normalization condition is, 



|dr(|G„,^(r)p + |F„,^(r)p) = l. 



(39) 



As usual, n specifies the angular quantum numbers and itp) the eigenvalue of the isospin 
operator (/|^). The total energy of the hypernuclear system is then given by 

+ 9u>^{r)pB{r) + ^5(p6(r)p3(r) + eA{r)pp{r)] 



- - I dr[-glCy{a{r))a{r)pl{r) 



+ 9l^{r)pl{r) + g],h{r)tppl{r) + eQy A{r) pl{r)] 



Y I 



(40) 



3 Spin-orbit potential in the QMC model 

In this section, we will focus on the spin-orbit potential of a hyperon in a hypernucleus. In 
order to see the differences between the QMC model and the QHD type models, we first discuss 
the spin-orbit potential of the A. 

The origin of the spin orbit force for a composite nucleon moving through scalar and 
vector fields which vary with position was explained in detail in Ref. [0] a) - cf. Section 3.2. 
The situation for the A is different in that, in an SU(6) quark model, the u and d quarks are 
coupled to spin zero, so that the spin of the A is carried by the s quark. As the a-meson is 
viewed here as a convenient parametrization of two-pion-exchange and the uj and p are non- 
strange, it seems reasonable to assume that the a, uj and p mesons couple only to the u and 
d quarks. The direct contributions to the spin-orbit interaction from these mesons (derived 
in Section 3 of Ref. a)) then vanish due to the flavor-spin structure. Thus, the spin-orbit 
interaction, V^q {r)l-s, at the position r of the A in a hypernucleus arises entirely from Thomas 
precession: 

Vio.(r)T- -Iv, X f . {±lMUr} + 9>(r)l) I- 1 (41) 
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where, v\ = p\/Mjl, is the velocity of the A in the rest frame of the A hypernucleus, and the 
acceleration, dv\/dt, is obtained from the Hamilton equations of motion applied to the leading 
order Hamiltonian of the QMC model . Because the contributions from the effective mass of 
the A, My^(r), and the vector potential, g^uj{r)^ are approximately equal and opposite in sign, 
we quite naturally expect a very small spin-orbit interaction for the A in the hypernucleus. 
Although the spin-orbit splittings for the nucleon calculated in QMC are already somewhat 
smaller [||, than those calculated in QHD p , we can expect much smaller spin-orbit splittings 
for the A in QMC, which will be confirmed numerically later. In order to include the spin- 
orbit potential of Eq. ( p]) correctly, we added perturbatively the correction due to the vector 
potential, — 2M*^{r)r ^ ' to the single-particle energies obtained with the Dirac 

equation Eq. (^, by evaluating it with the obtained wave function for the A. This is necessary 
because the Dirac equation corresponding to Eq. leads to a spin-orbit force which does not 
correspond to the underlying quark model, namely: 

This correction to the spin-orbit force, which appears naturally in the QMC model, may also 
be modelled at the hadronic level of the Dirac equation by adding a tensor interaction. 

In the QMC model, the general expression for the spin-orbit potential felt by the nucleon 
or hyperon, j (j = A, S, S), may be expressed as 



VL{r)l-s 
with 

Ai = 

G' = 



-1 



IJ'S 



2Mf{r)r I 



f-s, (43) 



d_ 

dr 



d ^1, 
dr 3' 



i=u,d 

< j\T.i=u4\^{i)\j > 



--<3\ E lfS)\3 > 



F] 



i=u,d 



-MnIo 



5 



IoMn 



3 2ilN{n]^ — 1) + rriu^R 



R] 



2m. 



u,d 



N 



:^)Rj - 3 



3 2fi^(a)(m(a)-l)-t-m:»i?*- 



(44) 
(45) 
(46) 

(47) 
(48) 

(49) 



Here r is the position of the baryon j in the hypernucleus (nucleus), and the terms proportional 
to /is and are the anomalous contributions from the u and/or d quarks due to the finite size 
of the hyperon (nucleon). They are related to the magnetic moments of the proton, /ip, and the 
neutron, /i„, as Us = l^p + f^n and /i^, = fip — fin, with the experimental values, fip = 2.79 and 
fin = —1.91 (in nuclear magnetons). However, these fig and /i^,, together with the quantities 
Iq and /* of Eqs. ( ^8] ) and (^9]), must be, and will be calculated self-consistently in the model. 
Note that the in-medium bag radius, R*, and the lowest bag eigenenergy for the u (d) quark. 
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n* in units of depend on j (j = A^, A, S, H), because the bag radius is generally different 

for each iso-multiplet of the octet baryons in the medium. The explicit expressions for the 
spin-orbit potentials for the octet baryons in the QMC model are given by: 



Vio.ir 



Vi 



s.o. 



-1 



-1 

2M*^{r)r 

-1 
2Mf{r)r 

-1 
2M^{r)r 

-1 
2M^{r)r 

-1 
2M^{r)r 

-1 
2Mf{r)r 

-1 



+ 3(1 - 2/i,r7^(r))A^ + ^(1 - 2^,r/^(r))A 



N 



1 



A;:'+3(l-2/i,r/^(r))A 
A:^ + 2A^ 



N 



1 - 2/i^?7iv(r))A^ 



2Mf( 



r r 



A^ + 2(1 - 4/x,r/s(r))AS + (1 - -^^,^^{r))^, 



A^ + 2(l-4/x,r/s(r))A^ 



A^ + 2(1 - 4^,772(r))AS - (1 - -/i,r7s(r))A, 
A= + (1 + 2/i,r7s(r))A= + ^(1 + ^^,r/H(r))A, 
Af + (1 + 2/i,r7H(r))A^ - ^(1 + '^^i,r]^{r))\ 



(50) 
(51) 
(52) 
(53) 
(54) 
(55) 
(56) 
(57) 



In the next Section, we will show the spin-orbit potentials for A, S'' and "E^ in yCa 
and y^Pb hypernuclei, which are evaluated self-consistently in the QMC model. These will 
be shown for the lp3/2 hyperon configuration, for which the change in the self-consistently 
calculated meson mean fields is expected to be the largest among the states which have non- 
zero spin-orbit interactions. 



4 Numerical results within the QMC model alone 

In this section, we will discuss some of the results calculated in the QMC model alone, in order 
to extract the modifications necessary to achieve a more realistic calculation. 

First, we need to specify the parameters used in the calculation. The parameters at 
the hadronic level, which are already fixed by the study of infinite nuclear matter and finite 
nuclei f^, are shown in Table ^. Concerning the parameters for the a meson, we note that the 



Table 3: Parameters at the hadronic level (for finite nuclei) 0. 



field 


mass (MeV) 


^V47r (eV47r) 


a 


418 


3.12 


ijj 


783 


5.31 


P 


770 


6.93 


A 





1/137.036 



properties of nuclear matter only fix the ratio, {ga/^a)-, with a chosen value, iria- = 550 MeV. 
Keeping this ratio to be a constant, the value = 418 MeV for finite nuclei (hypernuclei) is 
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obtained by fitting the r.m.s. charge radius of ^°Ca to the experimental value, rch(^''Ca) = 3.48 
fm||. 

Next, we show the calculated single-particle energies for y ^Pb hypernuclei in Table ^ In 



Table 4: Single-particle energies (in MeV) for y ^Pb hypernuclei for F = A, S and S. j^^Ph* 
denotes the results calculated with using the scaled coupling constant, 0.93 x g^{a = 0), 
which reproduces the empirical single-particle energy for the lsi/2 in ^^Ca, -20.0 MeV [|^. 
The experimental data for ^°^Pb are taken from Ref. ||3^. Spin-orbit splittings are not well 
determined by the experiments. 





f^Pb (Expt.) 


209 pi^* 




^^Ph 


|!i Pb 


zuy L) |-« 


§"Ph 


|i Pb 




-27.0 


-27.4 


-35.9 


-46.4 


-35.8 


-25.8 


-33.5 


-23.6 


lP3/2 




-23.8 


-31.9 


-41.7 


-32.0 


-22.6 


-29.7 


-20.5 


lPl/2 


-22.0 (Ip) 


-23.8 


-31.9 


-41.6 


-31.9 


-22.5 


-29.7 


-20.5 


14/2 




-19.5 


-27.1 


-36.8 


-27.6 


-18.6 


-25.6 


-16.8 


2Si/2 




-17.9 


-25.4 


-36.2 


-25.9 


-15.9 


-25.1 


-15.7 


lc?3/2 


-17.0 (Id) 


-19.5 


-27.1 


-36.6 


-27.4 


-18.3 


-25.6 


-16.8 


1/7/2 




-14.7 


-21.8 


-31.5 


-22.6 


-13.9 


-21.2 


-12.7 


2p3/2 




-12.6 


-19.4 


-29.7 


-20.3 


-10.9 


-20.4 


-11.3 


1/5/2 


-12.0 (1/) 


-14.6 


-21.7 


-31.1 


-22.3 


-13.5 


-21.3 


-12.8 


2pi/2 




-12.6 


-19.4 


-29.5 


-20.1 


-10.7 


-20.4 


-11.4 


lfl'9/2 




-9.5 


-16.0 


-25.8 


-17.2 


-8.6 


-16.7 


-8.3 


lfl'7/2 


-7.0 (Ig) 


-9.4 


-15.9 


-25.3 


-16.6 


-8.1 


-16.8 


-8.4 


l/ill/2 




-4.0 


-9.8 


-19.9 


-11.5 


-3.0 


-12.1 


-3.7 


2(i5/2 




-7.1 


-13.2 


-23.6 


-14.4 


-5.2 


-15.8 


-6.7 


2(^3/2 




-7.1 


-13.2 


-23.4 


-14.1 


-4.9 


-15.9 


-6.7 


1^9/2 




-3.9 


-9.7 


-19.2 


-10.7 


-2.2 


-12.2 


-3.8 


3Si/2 




-6.2 


-12.1 


-23.5 


-13.2 


-3.1 


-15.7 


-5.9 


2/7/2 




-1.8 


-6.9 


-17.7 


-8.4 




-11.6 


-2.0 


3p3/2 




-1.1 


-5.6 


-16.9 


-7.0 




-7.4 


-1.4 


2/5/2 




-1.7 


-6.8 


-17.3 


-8.0 




-11.7 


-2.0 


3pi/2 




-1.1 


-5.6 


-16.8 


-6.8 




-7.4 


-1.4 


1^13/2 






-3.4 


-13.9 


-5.4 




-7.6 





principle, the existence of the hyperon outside of the nuclear core breaks spherical symmetry, 
and one should include this in a truly rigorous treatment. We neglected this effect, since it 
is expected to be of little importance for spectroscopic calculations [^, |6^. However, we 
included the response of the nuclear core arising from the self-consistent calculation, which 
is significant for a description of the baryon currents and magnetic moments, and a purely 
relativistic effect [^, Thus, we should always specify the state of the hyperon in which the 
calculation was performed. 

Concerning the single-particle energy levels for A hypernuclei, although a direct compar- 
ison with the data is not precise due to the different configurations from those observed in the 
experiments, the calculated results seem to significantly overestimate the binding. In order 
to make an estimate of the difference due to this different configuration, we calculated the 
following quantity. By removing one IP3/2 neutron in ^^O, and putting a A as experimentally 
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observed, we calculate in the same way as for j^O - which means the nuclear core is still treated 
as spherical, although the core nucleus is deformed, and a rearrangement of the shell structure 
is expected to happen. In this case the calculated energy for the lsi/2 A is -19.9 MeV, to be 
compared with the value -20.5 MeV of the present treatment. 

At first look, the spin-orbit splittings for the A single-particle energies for ^°^Pb almost 
vanish, as explained in Section This is a very successful consequence of the explicit quark 
structure of the A considered in the SU(6) quark model. 

We repeated the calculation using the scaled coupling constant, 0.93 x g^{a = 0), which 
reproduces the empirical single-particle energy for the lsi/2 in ^^Ca, -20.0 MeV |Q. The results 
obtained using this scaled coupling constant are denoted by |°^Pb* in Table H. Then, one can 
easily understand that the QMC model does not require a large SU(3) breaking effect (only 
7 %) to reproduce the empirical single-particle energies. We should mention that a similar 
calculation was carried out by varying the coupling constant, g^^. In this case, the coupling 
constant obtained was, 1.10 x g^^, to reproduce the empirical single-particle energy for the lsi/2 
in j^Ca, -20.0 MeV The calculated results with this scaled coupling constant, 1.10 x gf^ 



A 



are almost identical to those calculated with the coupling constant, 0.93 x g^. This shows 
again that the QMC model does not require large SU(3) symmetry breaking to reproduce the 
empirical data. 

Although the present treatment is based on the underlying quark structure of the nucleon 
and hyperon, the overestimates discussed above may be ascribed to the lack of an explicit 
inclusion of the Pauli principle at the quark level among the u and d quarks in the core nucleons 
and the A. If one wants to reproduce the experimental single-particle energies, one needs to 
include the Pauli principle at the quark level which gives some repulsion. 

For S and S hypernuclei, the situation is even more complicated. We need to take into 
account the effect of channel coupling, SA^ — AA^ and EN — AA, as well as the Pauli blocking 
effect at the quark level. In particular, the Pauli blocking effect is significant in the open, 
coupled AN channel [|T5], |2|, |0|, |3|, ||, |6|, g^' H' 13 HI' whereas there are no such 
open channels for two nucleons in the core of a hypernucleus. The results shown in Table ^ 
(and also Table H) are obtained without including either these channel coupling effects, or the 
Pauli blocking effect at the quark level. These effects will be considered in specific ways at the 
hadronic level in Section 5. 

In Table 1^, we list the calculated binding energy per baryon (including hyperon), —E/A, 
r.m.s. charge radius, Vch, and r.m.s. radii of the the hyperon Y and the neutron and proton 
distributions (ry, r„ and r^, respectively), for the lsi/2 and 1^3/2 hyperon configurations. The 
r.m.s. charge radius is calculated by convolution with a proton form factor For comparison, 
we also give these quantities without a hyperon - i.e., for normal finite nuclei. The differences 
in values for finite nuclei and hypernuclei listed in Table ^ reflect the effects of the hyperon, 
through the self-consistency procedure. One can easily see that the effects of the hyperon 
become weaker as the atomic number becomes larger, and as the hyperon binding energies 
become smaller. We should comment on the small change in the r.m.s. radius of the protons, 
Tp, in Table |[ This appears to be a very subtle effect, related to the fact that the proton 
levels in the core are relatively shallow and closer to each other compared with those calculated 
in QHD 0. Thus, the proton wave functions in the QMC model are already pushed outside 
compared to those of QHD. We therefore expect then to be much more sensitive to the change 
in the mean field potentials and especially to the relative change of the Coulomb potential. We 
believe that this is the reason why Vp increases. For example, when the Coulomb potential in 
the core of the hypernucleus was artificially switched off (but the potential from the p meson 
was kept), both r„ and decreased for the lowest, lsl/2, hypernucleus state. 
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Table 5: Binding energy per baryon, —E/A (in MeV), r.m.s. charge radius, Teh-, and r.m.s. radii 
of the hyperon, ry, neutron, r„, and proton, (in fm) for y Ca and y^Pb. y^^Ca* and ^°^Pb* 
are the results calculated with using the scaled coupling constant, 0.93 x (7^(0" = 0), which 
reproduces the empirical single-particle energy for the lsi/2 in ^^Ca, -20.0 MeV |Q. Double 
asterisks, **, indicates the value used for fitting. 



hyper nuclei 


hyperon state 


-E/A 


Teh 


ry 






i^Ca* 




7.66 


3.51 


2.76 


3.31 


3.42 


A 


lSl/2 


7.83 


3.52 


2.55 


3.30 


3.42 


^+ Ca 


l5l/2 


7.65 


3.52 


2.69 


3.31 


3.42 


IjoCa 


lSl/2 


7.83 


3.52 


2.52 


3.30 


3.42 


IjLCa 


lSl/2 


7.97 


3.52 


2.31 


3.30 


3.42 


^^oCa 


lSl/2 


7.43 


3.50 


2.80 


3.32 


3.40 


^_Ca 


lSl/2 


7.54 


3.50 


2.62 


3.32 


3.40 


i^Ca^ 


lP3/2 


7.48 


3.50 


3.43 


3.32 


3.41 


A Ca 


lP3/2 


7.57 


3.51 


3.17 


3.32 


3.41 


^+ Ca 


lP3/2 


7.47 


3.51 


3.24 


3.32 


3.41 


^oCa 


lP3/2 


7.59 


3.51 


3.13 


3.32 


3.41 


E-Ca 


lP3/2 


7.70 


3.51 


3.02 


3.31 


3.41 


:^oCa 


lP3/2 


7.30 


3.49 


3.58 


3.33 


3.39 


HCa 


lP3/2 


7.40 


3.49 


3.36 


3.32 


3.39 


40Ca 




7.36 


3.48** 




3.33 


3.38 


209pb. 


lSl/2 


7.36 


5.49 


3.97 


5.67 


5.43 


^Pb 


lSl/2 


7.39 


5.50 


3.79 


5.67 


5.43 


|09pb 


lSl/2 


7.37 


5.49 


4.06 


5.67 


5.43 


|°o^Pb 


lSl/2 


7.39 


5.50 


3.77 


5.67 


5.43 


|09pb 


lSl/2 


7.41 


5.50 


3.42 


5.67 


5.43 


|?9pb 


lSl/2 


7.32 


5.49 


3.82 


5.68 


5.43 


2b9pj^ 


lSl/2 


7.34 


5.49 


3.48 


5.68 


5.43 


209pb^ 


lP3/2 


7.34 


5.49 


4.70 


5.68 


5.43 


209pb 


lP3/2 


7.37 


5.49 


4.50 


5.67 


5.43 


|09pb 


lP3/2 


7.34 


5.49 


4.65 


5.68 


5.43 


|o?Pb 


lP3/2 


7.37 


5.49 


4.46 


5.67 


5.43 


|09pb 


lP3/2 


7.39 


5.49 


4.28 


5.67 


5.43 


|Oo'Pb 


lP3/2 


7.30 


5.49 


4.57 


5.68 


5.43 


209pb 


lP3/2 


7.32 


5.49 


4.35 


5.68 


5.43 


208pb 




7.25 


5.49 




5.68 


5.43 
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Next we show in Figs. ^ and ^the spin-orbit potentials for the hyperons, A, S° and in 
y Ca and fP^Pb hypernuclei for the lp3/2 hyperon configuration. For comparison, we also show 
the spin-orbit potential for the neutron calculated self-consistently by Eq. (|5TD , with the same 
configuration as the hyperon. 

It is interesting to focus on some particular single-particle energy levels, e.g., lgg/2 and 
lg7/2 in |°^Pb, I'cPPb and lo^Pb. Comparing these single-particle energies with the plotted spin- 
orbit potentials in Fig. ^, one can understand the origin of the reverse in magnitudes for the 
single-particle energies between the I'cPPb (^°^Pb) and |o^Pb, and why the spin-orbit splittings 
for I'o^Pb are the most pronounced ones among them. Although the spin-orbit potentials shown 
in Fig. § are for the lp3/2 hyperon configuration, the main features also hold for these energy 
levels. 

As for the spin-orbit potentials for the S in hypernuclei, the QMC results show a slight 
inequality jVf'o | < jV^^ |, which is indicated in Fig. ^. This contrasts with the result of 
Pirner |jl3[ , | Vg'Q | = 1 1 Kso I ' which did not include the contribution from the Thomas preces- 
sion {Gj term in Eq. (^3])), but is rather similar to the results of Brockmann and Dover, 



Millener and Gal p6[. (Note that the ratio | of Ref. |]T3[ does agree with the ratio F^/F^ 
alone.) Finally, we should also note that the vector potential adopted by Pirner is based 
on the combined quark-gluon exchange. Thus, the corresponding spin-orbit potential of ours, 
F? term in Eq. (^3]), has a different origin. 



5 Effects of Pauli blocking and channel coupling 

As discussed in Section ^, when we insert a hyperon into an orbital with the same quantum 
numbers as one already occupied by nucleons we need to include the effects of Pauli blocking 
at the quark level. In the case of the E-hyperon there is an additional correction associated 
with the coupling to the open AN channel. These effects will be considered in specific ways at 
the hadronic level in this Section. A more consistent treatment at the quark level is beyond 
the scope of the present study. 

5.1 The Pauli blocking effect 

Here we consider the effect of the Pauli blocking at the quark level in a specific way at the 
hadronic level. This effect is significant in the open, coupled channel, particularly in the 
EA^ — AA^ channel in the present case, where such phenomenon simply does not occur for two 
nucleons in the core nucleus. Thus, the Pauli blocking effect we include in the present study 
should be regarded in order to ensure the correct channel coupling effect. 

It seems natural to assume that this effect works repulsively in a way that the strength 
is proportional to the u and/or d quark baryonic (number) density of the core nucleons. Then 
one can expect that the u and/or d quarks in the hyperon feel a stronger repulsion at the 
position where the baryon density is large. As a consequence, the wave function of the hyperon 
(quark) will be suppressed in this region. We assume here that the Pauli blocking effect is 
simply proportional to the baryonic density, although one could consider a more complicated 
density dependence. Then, the Dirac equation for the hyperon Y, Eq. (^), may be modified by 

[zj-d- M*{cr) - ( AypB(r-) + g^u{r) + (?p/|'6(f) + eQyA(f) ho]Mf^ = 0, (58) 

where, ps(r) is the baryonic density at the position r in the hypernucleus due to the core 
nucleons, and Ay is a constant to be determined empirically. (We note that the Pauli effect 
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associated with the hght quarks in the A should also lead to some repulsion for the nucleons. 
As this will be an 0{^) effect on the nucleon levels we have neglected this reaction back on the 
nucleons in the present work.) In the present treatment, we chose this constant Ay,y = A, in 
order to reproduce the empirical single particle energy for the lsi/2 in f'^Pb, -27.0 MeV [^. 



One might also imagine that this fitted value includes the attractive AA^ T^N channel 
coupling effect for the A single-particle enegies, because the value fitted is the experimentaly 
observed one. However, for E hypernuclei, the repulsive AN channel coupling effect 

must be included in addition to this effective Pauli blocking, in a way to reproduce the relative 
repulsive energy shift in the single-particle energies for the S. This issue will be discussed in 
the next Section. The fitted value for the constant Aa is Aa = 60.25 MeV (fm)^. Then for the 
E and E hypernuclei, we can deduce the constants, As,s, corresponding to the effective Pauli 
blocking effect as. As = Aa, and Ah = |Aa, by counting the total number of u and d quarks in 
those hyperons. 

5.2 The effect of channel coupHng 

We consider here the channel coupling (strong conversion) effect additional to the Pauli blocking 
at the quark level. It is expected that the channel couplings, SA^ — AA^ and SA^ — AA, generally 
exist in hypernuclei, and the former coupling is considered to be especially important [|l^, ^ 

H, 0, 13, |8|, g 



First, we consider the T,N — AN channel coupling. We estimate this effect using the 
Nijmegen potential |7I[] as follows. Including the effective Pauli blocking potential, AsPB(r) 
(As = Aa = 60.25 MeV (fm)^), we obtain as the single-particle energy for the lsi/2 level in 
209Pb, -26.9 MeV. This value does not contain the effect of the correct channel coupling to 
the A. On the other hand, in a conventional first-order Brueckner calculation based on the 



standard choice of the single-particle potentials (cf. ref. |]72[ for details) the binding energy for 
the S in nuclear matter with the Nijmegen potential fTlI] is, 12.6 MeV, for the case without 
including the correct channel coupling effect to the A, namely, without the Pauli-projector in 
the AA^ channel. When the channel coupling to the AA^ is recovered and the Pauli-projector 
in the AA^ channel is included, the binding energy for the S in nuclear matter decreases to 5.3 
MeV. Then, the decrease in the calculated binding energy for the S, 12.6 - 5.3 = + 7.3 MeV, 
should be considered as a net effect of the T^N — AN channel coupling for the S. In the present 
study, we include the effect by assuming the same form as that was applied for the effective 
Pauli blocking via AsP_B(r), and adjust the parameter As = Aa to As 7^ As to reproduce this 
difference in the single-particle energy for the lsi/2 in IjgPb, namely, -19.6 = -26.9 + 7.3 MeV. 
Here we should point out that our previous study 0] shows that the baryon density around 
the center of the ^osp]^ nucleus is consistently close to that of nuclear matter within the model. 
The value obtained for As in this way is. As = 110.6 MeV (fm)^. 

As for the EN — AA channel coupling, a study of Carr, Afnan and Gibson shows that 
the effect is very small for the calculated binding energy for %/Jie. The difference in the binding 
energy is typically less than 1 MeV for the results calculated with and without inclusion of the 
channel coupling effect. Although their estimate is not for large atomic number hypernuclei, 
nor nuclear matter, we expect that the EN — AA channel coupling effect is not so large for the 
calculation of the single-particle energies, and thus, we neglect it in the present study. 
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5.3 Results 



Including all the effects, the Pauli blocking and the UN — AA^ channel coupling discussed in 



Sections |5.1| and |5.2| , we will present the final results in this Section. 

In Tables |^ and 0, we list the calculated single-particle energies for O, y^Ca, y^Ca, y^Zr 
and y^Pb hypernuclei, together with the experimental data for the A hypernuclei. Unfortu- 
nately, data for the larger atomic number hypernuclei are limited to A hypernuclei. Concerning 
the single-particle energy levels for the A hypernuclei, the QMC model supplemented by the 
effective Pauli blocking effect employed at the hadronic level, reproduces the data reasonably 
well. The small spin-orbit splittings are still achieved, as before, with the QMC model alone. 

It is interesting to compare the single-particle energies for the charged hyperons, and 
S~, and those of the neutral hyperons, S'' and rP. The present results imply that the Coulomb 
force is important for forming (or unforming) a bound state of the hyperon in hypernuclei. This 
was also discussed by Yamazaki et al. [p5|, in the context of the light S~ hypernuclei. 



Table 6: Single-particle energies (in MeV) for y O, y Ca and y Ca hypernuclei, calculated with 
the effective Pauli blocking and the SA^ — AA^ channel coupling. Experimental data are taken 
from Ref. fS^. Spin-orbit splittings are not well determined by the experiments. 



O (Expt. 



A ^ 



17 

s- 



O 



17 o 



17 O 



17 



o no 



lSl/2 
lP3/2 
lPl/2 



-12.5 
-2.5 (Ip) 



-14.1 
-5.1 
-5.0 



-17.2 
-8.7 
-8.0 



-9.6 
-3.2 
-2.6 



-3.3 



-9.9 
-3.4 
-3.4 



-4.5 



i^Ca (Expt. 



Ccl 



^0 Ca 



41 p 41 p 41 p 



lSl/2 
lP3/2 
lPl/2 



Id. 
2s 



5/2 
1/2 

3/2 



-20.0 
-12.0 (Ip) 



-19.5 
-12.3 
-12.3 
-4.7 
-3.5 
-4.6 



-23.5 
-17.1 
-16.5 
-10.6 
-9.3 
-9.7 



-13.4 
-8.3 
-7.7 
-2.6 
-1.2 
-1.9 



-4.1 



-17.0 
-11.2 
-11.3 
-5.5 
-5.4 
-5.6 



-8.1 
-3.3 
-3.4 





^ Ccl 


f-Ca 


^0 Ca 


If+Ca 


l-Ca 


foCa 


lSl/2 


-21.0 


-19.3 


-14.6 


-11.5 


-14.7 


-12.0 


lP3/2 


-13.9 


-11.4 


-9.4 


-7.5 


-8.7 


-7.4 


lPl/2 


-13.8 


-10.9 


-8.9 


-7.0 


-8.8 


-7.4 


14/2 


-6.5 


-5.8 


-3.8 


-2.0 


-3.8 


-2.1 


2Si/2 


-5.4 


-6.7 


-2.6 




-4.6 


-1.1 


14/2 


-6.4 


-5.2 


-3.1 


-1.2 


-3.8 


-2.2 


1/7/2 




-1.2 











In Table ^, we show the calculated binding energy per baryon, r.m.s charge radius, r^h 
and r.m.s. radii of the hyperon, neutron and proton distributions, ry, r„ and r^, respectively. 
The results listed in Table |] are calculated with the lsi/2 hyperon configuration for all cases. 
Here, one can again notice the important role of the Coulomb force. For example, the binding 
energy per baryon, —E/A, for the S" (S+) hypernuclei is typically the largest (smallest) among 
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Table 7: Single-particle energies (in MeV) for y Zr and y^Pb hypernuclei, calculated with the 
effective Pauli blocking and the TjN — AN channel coupling. Experimental data are taken from 
5B| . Spin-orbit splittings are not well determined by the experiments. Double asterisks, 



Ref. 



indicates the value used for fitting. 





^^Yb (Expt.) 


91 


91 ry 


91 


91 

E+Zr 


91 '7 

|_Zr 


91 '7 

|{,Zr 




-22.5 


-23.9 


-27.3 


-16.8 


-8.1 


-22.7 


-13.3 


-1 

lP3/2 




-18.4 


-20.8 


-12.7 


-5.0 


'\ ''7 A 

-17.4 


-9.7 


-1 


-lo.U (Ipj 


-18.4 




-12.4 


A '1 

-4.7 


A 'T A 

-17.4 


-9.7 


1 ^ 




-iz.o 


1 /I 
-I0.4 


-o.i 


-U.y 


-Iz.o 


-0.4 


2Si/2 




-10.8 


-15.6 


-6.5 




-12.4 


-3.9 


14/2 


-9.0 (let) 


TOO 

-12.3 


-14.8 


-7.5 


-0.3 


-12.4 


-5.5 


1/7/2 




-5.9 


-10.2 


"1 

-0.1 
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the same atomic number hypernuclei, while the r.m.s. radii for the hyperon, ry, r^- (^"2+) is 
mostly the smallest (largest) among the hypernuclei of the same atomic number. 

Finally, we show the strengths of the scalar and vector potentials for the hyperons in y O, 
y Ca and y ^Pb with y = A, and H*^ in Figs. and the effective masses of the nucleon and 
hyperon together with the baryon densities in ^°^Pb, |*(?Pb and lo^Pb in Figs. ^ /^From the 
baryon densities given for each hypernucleus, one can estimate the magnitude of the effective 
Pauli blocking and the — AN channel coupling effect employed in the present study. Note 
that the baryon densities versus the distance r from the center of the hypernuclei are slightly 
different for each hypernucleus due to the self-consistent calculations. We should also mention 
that the results shown in Figs. ^ are calculated by employing different parameters (and a 
different version of the model) from those shown in Ref . . 



6 Summary and discussion 

In summary, we have reported results for A, S and S hypernuclei which were systematically 
calculated with the QMC model for the first time. Especially, the spin-orbit potentials for the 
hyperons in the hypernuclei were evaluated self-consistently with the explicit quark structure 
for the bound hyperon and nucleons. The very small spin-orbit force for the A in hypernuclei 
was achieved naturally in the present treatment. This is a direct consequence of the SU(6) 
quark model wave function for the A used in the QMC model. However, the single-particle 
energies calculated in the QMC model alone tend to overestimate the experimental data. In 
order to overcome this overbinding problem, we included the effect of Pauli blocking in a 
phenomenological way. For the S hypernuclei, we took into account the channel coupling 
effect estimated using the Nijmegen potential. In the future, these effects, which were included 
phenomenologically at the hadronic level in the present study, should be treated on the same 



footing, namely, at the quark level. For this purpose, it seems that the quark cluster model [19 



20| may offer some possibilities to include the Pauli blocking effect at the quark level. However, 
an application of the model to larger atomic number hypernuclei, or nuclear matter, seems to 
be difficult, and has not yet been carried out. 

Although, there seem to be several points which could be improved in the present treat- 
ment, we would like to emphasize that, this study is the first, systematic investigation of the 
properties of hypernuclei based on explicit quark degrees of freedom. One of the advantages 
of the present study is the relatively simple nature of the QMC model, and simultaneously, its 
capability to give a quantitative description of the properties of finite nuclei. 
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Table 8: Binding energy per baryon. —E/A (in MeV), r.m.s charge radius, Teh-, and r.m.s radii 
of the hyperon, ry, neutron, r„, and proton, (in fm) for y O, y Ca, y Ca, y Zr and y^Pb 
hypernuclei. The configuration of the hyperon, Y , is lsi/2 for all hypernuclei. For comparison, 
we list also the corresponding results for normal finite nuclei. Double asterisks, **, indicates 
the value used for fitting. 



hypernuclei —E/A Vch ry r„ Vp 

fO K37 2^84 2.49 2.59 2.72 

i7+0 5.91 2.82 3.15 2.61 2.70 

^"^oO 6.10 2.83 2.79 2.60 2.71 

^^_0 6.31 2.84 2.49 2.59 2.72 

go 5.86 2.80 2.98 2.62 2.68 

6.02 2.81 2.65 2.61 2.69 

^^O 5.84 2.79 — 2.64 2.67 

^^Ca 7.58 3.51 2.81 3.31 3.42 

IVCa 7.33 3.50 3.43 3.32 3.41 

fiCa 7.44 3.51 3.14 3.31 3.41 

f-Ca 7.60 3.51 2.87 3.31 3.41 

foCa 7.34 3.49 3.09 3.32 3.40 

|lCa 7.45 3.50 2.84 3.32 3.40 

^"Ca 7.36 3.48** — 3.33 3.38 

fCa 7^58 3^54 2.84 3.63 3.45 

f+Ca 6.32 3.57 3.47 3.71 3.48 

foCa 7.40 3.54 3.14 3.64 3.45 

f_Ca 7.48 3.55 2.60 3.63 3.45 

i^oCa 7.32 3.53 3.14 3.65 3.43 

f-Ca 7.35 3.53 2.79 3.65 3.44 

^^Ca 7.27 3.52 3.66 3.12 

fzi 7^95 429 3.25 4.29 4.21 

IVZr 7.82 4.28 4.01 4.30 4.20 

|^,Zr 7.87 4.29 3.56 4.30 4.21 

|lZr 7.92 4.29 2.89 4.29 4.21 

|{,Zr 7.83 4.28 3.54 4.30 4.20 

l-Zr 7.87 4.28 2.98 4.30 4.20 

^°Zr 7.79 4.27 — 4.31 4.19 

f^Fh 7.35 5.49 3.99 5.67 5.43 

|?Pb 7.28 5.49 4.64 5.68 5.43 

|Oo^Pb 7.31 5.49 4.26 5.67 5.43 

|l9pb 7.34 5.49 3.80 5.67 5.43 

fo^Pb 7.30 5.49 3.96 5.68 5.43 

f^Fh 7.32 5.49 3.59 5.68 5.43 

^bsp^ 7^25 K49 — 5.68 5.43 
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Figure 3: Spin-orbit potentials, V^f^ (r), calculated in the QMC model alone for Y 
and S° in y Ca hypernuclei for the lp3/2 hyperon state. 
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Figure 4: Spin-orbit potentials, V^f^ (r), calculated in the QMC model alone for Y 
and S° in fP ^Pb hypernuclei for the lp3/2 hyperon state. 
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Figure 5: Potential strengths, g^{cr) and , for F = A, S° and S'^ in yO, yCa and y^Pb 
hypernuclei for the lsi/2 hyperon state, calculated with including the effects of Pauli blocking 
and the Y:N — AN channel coupling. 
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Figure 6: Effective masses of the nucleoli and hyperon, and the baryon density in y ^Pb hyper- 
nuclei for Y = A,I1^ and S°. The hyperon configuration is the lsi/2 state for all cases. 
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